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Abstract 

When recent experimental positronium (Ps) formation cross sections in noble gases 
have been compared with the most up-to date theoretical studies, the agreement is 
qualitative, but not quantitative. In this paper we re-examine this process and show 
that at low energies Ps formation must be treated nonperturbatively. We also look 
at Ps formation with inner shell electrons. 
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1 Introduction 



Positronium (Ps) represents a bound state between a positron and an electron. 
It is formed in positron-atom collisions, 

A + e + — >A+ + Ps, (1) 



when the positron energy, e = k 2 /2, is above the Ps formation threshold, 
e > \e n \ - \E U \ (2) 



where e n is the energy of the bound electron atomic orbital n, E\ s »s —6.8 eV 
is the energy of the ground-state Ps, and k is the incident positron momentum 
(atomic units are used throughout). 
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Recently positronium formation in Ne, Ar, Kr and Xe has been determined 
by two experimental groups [1,2]. The two sets of data are in fairly good 
agreement, especially at lower energies. However, recent distorted- wave Born 
approximation (DWBA) calculations [3] overestimate the cross sections by a 
large factor, ranging from 1.6 in Ne to 3 in Xe, while the overall energy depen- 
dence of the DWBA cross sections is in reasonable accord with experiment. 
This is in contrast with earlier coupled-static calculations [4], which yield 
better magnitudes of the cross section maxima, but disagree on the energy 
dependence. 

In this paper we perform lst-order and all-order calculations of Ps formation 
from valence and subvalence subshells. Our consideration is restricted to Ps 
formation in the ground-state. Noble gas atoms have tightly bound electrons, 
making excited-state Ps formation much less probable (see, e.g., [3]). We argue 
that a structure observed at energies beyond the main cross section maximum 
(described as a shoulder, or in some cases seen as a secondary peak [1,5]) 
is most likely related to Ps formation by the subvalence ns electrons. We 
also consider Ps formation from inner shells. It produces inner-shell vacan- 
cies and can be important for positron-annihilation-induced Auger-electron 
spectroscopy [6]. 



2 lst-order approximation 

2. 1 Ps formation amplitude and cross section 

Using lst-order many-body perturbation theory, and neglecting the interaction 
between the outgoing Ps and residual ion, the Ps-formation amplitude can be 
written as [7] 

(^i s , K \V\n,e) = J *L, K (ri,r 2 ) f- -— L— - j ^2)^1)^1*2, (3) 



where (f £ is the incident positron wavefunction, ip n is the Hartree-Fock wave- 
function of the initial electron state ("hole"), and ^i Sj k is obtained from the 
wavefunction of the ground-state Ps with momentum K, 

*i aiK (ri, r 2 ) = e iK -( ri+r2 )/ 2 ls ( ri - r 2 ), (4) 
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by orthogonalising it to all electron orbitals n' occupied in the target ground 
state, 



The positron wavefunction is calculated in the field of the ground state atom 
described in the Hartree-Fock approximation. The Ps center-of-mass motion 
is described by a plane wave. The Ps formation cross section is found by 
integration over the directions of K, 



where M = 2 and K = [2M(e — \e n \ + |-Ei s | )] 1//2 are the Ps mass and momen- 
tum, and e (r) ~ e* k r normalization is assumed. The approximation (3)-(6) 
is equivalent to DWBA for a rearrangement collision. 

A numerical calculation of the amplitude and cross section was performed by 
expanding the Ps wavefunction in electron and positron spherical harmonics 
with respect to the nucleus. Integration over the angular variables was done 
analytically, while the radial integrals are calculated numerically (see [7] for 
some details). To ensure accurate positions of the Ps formation thresholds, ex- 
perimental ionization energies \e n \, rather than the Hartree-Fock values, were 
used in the calculations. Detailed below are the Ps-formation cross sections for 
neon, argon, krypton and xenon, calculated using the lst-order approximation 
described above. 



2.2 Partial-wave contributions 

The cross sections are found by summing over the positron partial waves from 
I = to 10. Figure 1 shows the partial wave contributions for the 3p subshell 
of argon. Note that the p, d, f and g waves have the largest individual cross 
sections and make up most of the cross section peak. The contributions of 
higher partial waves are suppressed by the centrifugal barrier, preventing the 
close encounters which lead to Ps formation. The small contribution of the 
s-wave is due to it being spherically symmetric, making it harder for the 
positron to bind and move away with an electron. This is true for all the 
noble gases. The exceptionally small s-wave contribution to Ps formation was 
noticed earlier for hydrogen and helium (see, e.g., [8,9,10]) and explained by 
the hidden- crossing method [11]. 




(5) 




(6) 
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Positron energy (eV) 

Fig. 1. Partial wave contributions to the total Ps formation cross section for the 3p 
subshell of argon. Various thin curves are the contributions of I = 0-10, while the 
solid thick curve is the total cross section. 

2.3 Comparison with experiment 

Figure 2 shows the Ps-formation cross sections for the valence np and subva- 
lence ns orbitals together with their sum, for Ne, Ar, Kr and Xe. The present 
results for the np subshell practically coincide with the Ps(ls) formation cross 
section from DWBA [3]. The calculations are compared with the experimental 
data obtained with a cold, trap-based positron beam [2], and with the cross 
section found by subtracting the direct ionization cross section from the total 
ionization cross section [1]. Moving from Ne to Xe, the calculations increas- 
ingly overestimate the measured cross section near the maximum. 

For Ne, Ar and Xe experiment and theory converge at higher energies, while 
in Kr the discrepancy persists. Such convergence should be expected from a 
theory point of view. Indeed, at higher positron energies the dominant contri- 
bution to the amplitude (3) and cross section (6) comes from higher partial 
waves, for which the plane-wave description of the Ps motion is more accu- 
rate. At the same time, the contributions of individual partial waves to the 
amplitude become small. This means that higher-order corrections neglected 
by the lst-order theory may not be important (see below). Thus, we cannot 
offer an explanation for the divergence between theory and experiment in Kr. 

2-4 Inner-shell Ps formation 

Ps-formation thresholds for the inner shells lie at much higher energies, e.g., 
at 242 and 320 eV for the 2p and 2s orbitals in Ar. As a result, the incident 
positron wavefunction oscillates rapidly, reducing the magnitude of the am- 
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Fig. 2. Comparison of the lst-order Ps-formation cross sections for Ne, Ar, Kr 
and Xe with experiment. Calculations: long-dashed curve, contribution of the np 
subshell; short-dashed curve, contribution of the ns subshell; solid curve, total cross 
section; dotted curve, DWBA [3] for the np subshell. Experiment: solid circles, 
University of California at San Diego (UCSD) [2]; open diamonds, University College 
London (UCL) [1]. 



plitude (3). Ps formation from inner shells is additionally suppressed by the 
positron repulsion from the nucleus. The Ps-formation cross sections for Ne, 
Ar, Kr and Xe are shown in Figure 3. It can be noted that significantly higher 
energies are required to produce positronium for the lighter noble gases (with 
the more tightly bound electrons). 

In general, the Ps-formation cross sections for the inner-shell electrons are 
quite small. Thus, in neon, the valence shell contribution still dominates, even 
at positron energies of 1600 eV. In argon, it can be seen that up to 500 eV, the 
valence contribution dominates the cross section. However, above this energy 
all subshells contribute approximately equally to the total cross section, and as 
the energy increases further, the 2p subshell contributes the most. In krypton, 
the 3d subshell contributes most to the Ps-formation cross section, as this 
subshell is relatively far away from the nucleus. This is also true for xenon, 
where the Ad subshell dominates the cross section. Note also that in Kr and 
especially Xe the inner-shell cross sections are much larger than those of Ne 
and even Ar, with their Ps formation threshold values nearly ten times smaller 
than that of Is in Ne. 
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Fig. 3 
Xe. 



Ps-formation cross sections for inner and valence subshells of Ne, Ar, Kr and 
3 Ps formation: nonperturbative approach 



3. 1 Check of unitarity 



As seen from Figure 2, the Ps formation cross sections increase dramatically 
from Ne to Xe. This increase is matched by a growing discrepancy between 
the lst-order results and experiment. This suggests that as the cross sec- 
tions become larger, the lowest-order perturbation theory treatment becomes 
increasingly inaccurate. Indeed, it turns out that in the lst-order approxima- 
tion, equations (3)-(6), the lower partial- wave contributions (I = 0-3) which 
dominate near the cross section maximum, become close to and even violate 
(for Kr and Xe) the unitarity limit for inelastic processes, a® < n(2l + l)/k 2 . 
A comparison with the unitarity limit is presented in Figure 4 for the s and p 
partial waves in Ar, Kr and Xe. 

Physically, unitarity ensures that the amount of Ps formed in positron-atom 
collisions cannot be greater than the number of positrons going in. Any in- 
elastic cross section above this limit is physically impossible. The fact that 
our lst-order results (and the analogous DWBA cross sections of Ref . [3]) are 
close to, or exceed the unitarity limit means that Ps formation is a strong 
process which cannot be treated perturbatively. Hence, in addition to the 1st- 
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Fig. 4. Comparison of the s- and p-wave partial Ps formation cross sections for Ar, 
Kr and Xe with the unitarity limit. Thin curves, lst-order approximation; thick 
curves, all-order approximation; dotted curves, unitarity limit, ir(2l + l)/k 2 . 

order Ps-formation amplitude (3), one must include higher-order contributions 
which account for the decay of A + Ps back into the e + + A channel, followed by 
Ps formation, etc. In other words, the effect of Ps formation on the incident 
positron must be taken into account. We do this by means of an all-order 
approach outlined in the next section. 



3.2 All-order approximation 



The effect of Ps formation on the positron can be described by the Ps- 
formation contribution to the positron-atom correlation potential, defined by 
its matrix elements [7], 

/^(Ps)u\_ f W,n\V\^ ls ^)(^ ls ^\V\n,e) al 3 K 
{£\^e \e)-J E + £n _ Eu _ K 2 /A + l0 (27r) 3> 



where {^is,K.\V\ n ,£) is the amplitude (3), E ls + K 2 /4 in the denominator is 
the Ps energy, and the integral is over all Ps momenta K. Note that in this 
section we use positron states with a given angular momentum, i.e., spherical 
waves, (p E {r) = r~ 1 Yi m (Q)P e i(r), with the radial wavefunctions normalized by 
P £ i(r) ~ (7r/c) _1//2 sin(/cr — ^7rZ + ^ HF ), where (5 ; HF is the positron phaseshift in 
the static field of the Hartree-Fock atom. In this case the correlation potential 
(7) is determined separately for each positron partial wave. 

Below the Ps-formation threshold (e'|Sg S ' ) |e) is real. Above the Ps-formation 
threshold, for E > \e n \ — \Ei 8 \, the correlation potential acquires an imagi- 
nary part. This gives rise to "absorption" of the positron flux, which is being 
redirected into the Ps formation channel. In fact, the lst-order Ps formation 
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cross section (6) is proportional to the imaginary part of (e|S( Ps )|e), 



Im 



(s\^\s) = -rr^ J \(^ ls , K \V\ne) 



(8) 



where K = 2{e — \e n \ + l-E^I) 1 / 2 . Because of the different normalization of 
the positron states adopted in Sees. 2 and 3.2, caution is required when using 
equations (6) and (8) to relate a Ps to Im(e|S( Ps ) \e) (such relation is derived 
below) . 

A nonperturbative ( "all-order" ) calculation of Ps formation is done by solving 
an integral equation for the matrix elements of E s (see, e.g., [12]), 

(e'\t E \e) = (£'|4 Ps) k> + / W^WVP™^ ^ ( 9 ) 

The positron scattering phaseshift is then obtained as 

8 t = 5[ + i8'{ = 5f F + A5i (10) 



with 



tan A5i 



-7r(e|I! e |e) 



where A5i is the additional phaseshift due to the correlation potential. For 
energies E above the Ps formation threshold, the correlation potential, £^ >s ' ) , 
is complex, and the phaseshift has a nonzero imaginary part, 6" > 0. The Ps- 
formation cross section is then obtained from 8" by summing over the partial 
waves, 

oo 

a Ps = ^E(2^+l)(l-e- 45 "). (12) 



If we assume that Ps formation is a weak process, then (e|S^ >s ' ) |5) is small, and 
we have (£|E E |£) m (e|S^ >s ' ) |£), A8 t m — n(e\fl e \e), and 8" — 7rIm(e|S( Ps )|e). 
Using 5" < 1, we then have from equation (12): 



a 



(0 

Ps 



47T 



2 



k 2 



■(2/ + l)Im(e|S( Ps )|e). (13) 



The right-hand side of this equation is equivalent to the lst-order approxima- 
tion examined earlier in this paper. 
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3.3 Total Ps-formation cross sections 



Figure 5 shows the all-order and lst-order cross sections for Ne, Ar, Kr and 
Xe along with the two sets of experimental data. As expected, the differ- 
ence between the all-order and lst-order results is greater for atoms where 
the Ps formation cross section is large, i.e., for Kr and Xe. The all-order ap- 
proximation has reduced the cross section maxima for all atoms, but still not 
significantly enough to match the experiment. However, theory and experi- 
ment are in better agreement at higher positron energies. The situation looks 
especially encouraging in Xe, where both experiments and theory are close 
above 40 eV. Note that for Kr and Xe, the all-order and lst-order cross sec- 
tions are markedly different even at the higher-energy part of the scale. This 
emphasizes the need for nonperturbative treatment of Ps formation. 
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Fig. 5. Comparison of the calculated all-order and lst-order Ps-formation cross sec- 
tions for Ne, Ar, Kr and Xe with experiment. Calculations: dotted curve, total 
(lst-order); short-dashed curve, ns subshell (1st order); long-dashed curve, np sub- 
shell (all-order); chain curve, ns subshell (all-order); solid curve, total cross section 
(all-order). Experiment: solid circles, UCSD [2]; open diamonds, UCL [1]. Vertical 
dashed lines indicate Ps-formation thresholds for the subvalence ns electrons. 

Regarding the contribution of the subvalence ns subshell, the corresponding 
Ps-formation thresholds for Ne, Ar, Kr and Xe are 41.67, 22.47, 20.71 and 
16.59 eV, respectively 1 . The all-order calculation was not performed for Ne 



Note that the ns thresholds indicated in figures 5-8 of Ref. [1] are ns ionization 
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and Ar, where we expected its effect to be insignificant. In fact, the 2s-subshell 
contribution to the Ps-formation cross section in Ne is very small. Experimen- 
tal data in Ne also do not reveal any clear features that could be related to 
the opening of a new channel at e ~ 42 eV. The Ps formation from the 3s 
subshell in Ar is relatively more important. The opening of this channel at 
about 22.5 eV coincides with the onset of the second peak in the UCL data 
[1]. It also marks the start of a weak shoulder-like structure in the UCSD data 
[2], where earlier experiments [5] showed a much more prominent feature. The 
reason for the discrepancy between different experimental observations is at 
present unclear, but we suggest that the origin of this structure is most likely 
related to the opening of the 3s Ps-formation channel. 

Using the all-order approach makes a much greater difference in krypton and 
xenon. The valence np contribution is greatly reduced in the energy range 
around the maximum. In contrast, the subvalence ns contribution is increased 
in the all-order approach. The onset of the Ps-formation cross section from 
the 4s orbital in Kr is very smooth. The UCSD data for Kr do not show any 
feature just above 20 eV, while the UCL data possess a clear shoulder, whose 
onset is more rapid than that predicted by the theory. At higher energies 
the two sets of experimental data diverge, with our all-order results being in- 
between. Given that the theory is expected to be more reliable here than at 
lower energies, we suggest that this discrepancy points to a need for further 
experimental studies. 

Of the four atoms examined, Xe has the largest Ps formation cross section 
by the subvalence electrons. Its contribution results in a change of slope of 
the calculated total cross section. For Xe both sets of experimental data show 
a clear shoulder-like structure, whose onset is close to the 5s Ps-formation 
threshold. The overall size of the shoulder is comparable with the calculated 
5s Ps formation cross section (chain curve in Figure 5). At higher energies 
there is good agreement between the two experiments and the calculated cross 
section. 

One may speculate that in a better calculation, the Ps-formation cross section 
from the valence np orbital will be suppressed around its maximum and at 
the energies below 40 eV. Then, even if the ns Ps-formation cross sections 
remain close to the present estimates, their contribution to the total will be 
more noticeable. 

Finally, an alternative explanation of the shoulder/secondary-peak structure 
discussed, e.g., in Ref. [1], is that it is caused by Ps formation in excited states. 
According to the DWBA calculations of Gilmore et al. [3], this contribution 
is not negligible (though small). However, the thresholds for excited-state Ps 
formation lie much lower than the energies where the structures are observed, 

thresholds, rather than the Ps-formation thresholds. 
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making its importance for these structures questionable. 



4 Conclusions 



A comparison of the lst-order and all-order results shows that as Ps forma- 
tion is strong, it cannot be treated perturbatively. Going beyond the lst-order 
approximation reduces the cross sections, especially at low energies. However, 
below 40 eV the calculated cross sections are still higher than experimen- 
tal values. Above this energy theory and experiment generally converge. In 
particular, in Xe we observe good agreement between the calculations and 
experimental data from the UCL and UCSD groups. 

There are two reasons for the discrepancy between theory and experiment 
at low energies. First, the motion of Ps in our calculations is described by 
a plane wave. The electron part of the Ps wavefunction is orthogonalized 
to the target electron orbitals. This manifestation of the Pauli principle to 
some extent accounts for the interaction between the Ps and the final-state 
ion. On the other hand, the positron repulsion from the nucleus in the final 
Ps state is completely neglected. This repulsion is especially important for 
the lower positron partial waves. Its neglect is probably the main reason for 
the overestimation of the Ps-formation cross section maxima by the present 
method. 

Secondly, all open channels, i.e., elastic scattering, Ps formation and direct ion- 
ization must be included simultaneously. Above the atomic ionization thresh- 
old all of these channels compete for the positron flux. This effect can be 
accounted for by the correlation potential method described in Sec. 3.2, by 
adding the lowest 2nd-order contribution to (see, e.g., [13]). However, to 
be able to extract the Ps formation cross section from such a calculation, the 
formalism of Sec. 3.2 must be extended. In the present form the imaginary 
part of the phaseshift allows one to find only the total reaction cross section. 

One may expect that both of these effects will make the Ps formation cross 
section smaller, and bring it into a close agreement with experiment. By further 
including Ps formation from the inner valence subshell and Ps formation in 
excited states, one should achieve a complete description of the Ps formation 
process, including any secondary structures. 
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